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Abstract. A new parametrization (one-to-one onto map) of compact wavelet matrices of 
rank m and of order and degree N is proposed in terms of coordinates in the Euclidian space 
(j-<(r?i-i)Af _ rpj^g developed method depends on Wiener-Hopf factorization of corresponding 
unitary matrix functions and allows to construct compact wavelet matrices efficiently. Some 
applications of the proposed method are discussed. 
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1. Introduction 

Let C be the set of complex numbers, and let F C C be a subfield invariant under 
the complex conjugation. A formally infinite matrix A with m rows 



A 



(■■■ 



m— 1 
■^0 



1 



7 



, a) G F, 



is called a wavelet matrix (of rank m) if its rows satisfy the so called shifted orthogo- 
nality condition: 



(2) 



{5 stands for the Kronecker delta). Such matrices are a generalization of ordinary 
mxm unitary matrices and they play the crucial role in the theory of wavelets [7] and 
multirate filter banks [8]. The reason for us to work with F instead of just C is that 
it will allow more flexibility in applications, since the proposed proofs and discussions 
apply to a whole range of fields including the set of real numbers, rational numbers, as 
well as algebraic extensions of rational numbers. A reader not concerned with general 
fields may assume that F = C. 

In the polyphase representation of matrix A, 



(3) 



A(z) 



where A 
Ak = (a 



■ ■ ■ , v4_i, Ao, Ai, A2, ■ ■ ■) is the partition of A into block mxm matrices 
fcm+j)' 0<z,j<m — 1, the condition ([2]) is equivalent to 



(4) 



A(z)A(z) 
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where A{z) = Yl'k'=-oo ^fc-^"'^ adjoint to A{z) {A* := A is the Hermitian 

conjugate and Im stands for the m x m unit matrix). This is easy to see as (|2]) can be 
written in the block matrix form YlT=-oo ^kA^_f. = Siolm- 

Our notion of a wavelet matrix is weaker than usual. So, as the orthogonal basis of 
L^(M) can be developed (see [71 Ch-s 4, 5]), also the linear condition A(l) 1 = i/mei, 
where 1 = (1, 1, ... , 1)-^ and ei = (1, 0, . . . , 0), must be satisfied. In our consideration, 
the linear condition is irrelevant. Instead, we require the condition 

(5) A(l) = Im 

which much simplifies the whole exposition. For any wavelet matrix A, A(l) = 
'^kL-oo^k is unitary (see (jl])), so that the conditions ([2]) and ([5]) will be satisfied 

for ^0 = (^(1)) Note also that for any wavelet matrix satisfying (jS]) and any 
unitary matrix U, the matrix A = UAq is also a wavelet matrix satisfying A(l) = U. 
Thus the additional restriction does not lose generality in the description of wavelet 
matrices. Observe that the polyphase representation of UA is YlT=-oo UA^z'' (obvi- 
ously, there is one-to-one correspondence between the matrices ([1]) and their polyphase 
representations ([3]) and they are naturally identified). 

We consider the compact wavelet matrices, which means that only a finite number 
of coefficients in ([T]) are non-zero (the corresponding wavelet functions in L^(]R) have 
then a compact support, and the corresponding filters in signal processing applica- 
tions are physically realizable). Namely, compact wavelet matrices with polyphase 
representation 

N 

(6) A{z) = Y,A,z\ 

k=0 

where An 7^ 0, are called of (rank m and) order N := ord(^) (in some works they are 
called of genus + 1); we write A G WA^(m, A^, F). The property (jl]) for the matrix 
polynomial (|6]) means that A{z) is a paraunitary matrix function. Note that, in this 
case, A(z) is usual unitary matrix for each z eT := {z E : \z\ = 1}. 

As the determinant of A{z) G VVA^(m, A^, F) is always monomial, i.e., det A(z) = 
cz'^ (since det A(2;) det A(z) = 1), the positive integer d is defined to be the degree 
oi A, d = deg(^). The set of wavelet matrices of rank m, order N and degree d 
will be denoted by >VA^(m, A^, ci, F). It can be proved that deg(^) > ord(^) for 
any A G VVA^(m, A^, F) and deg(^) = ord(^) if and only if rank(Ao) = m — 1 (see 
Lemma 1 below). As AqA*^ = for each A G WA^(m, A^, F) (since it is the A^th 
matrix coefficient of the product in (^) and 7^ 0, we have rank(74o) < m. Thus 
deg(^) = A^ whenever Aq has a maximal possible rank and it is assumed to be the 
nonsingular case. Throughout the paper, we will intensively study WA^(m, A^, A^, F). 

Obviously, if we shift rows of ^ G WA^(m, A^, A^, F) by arbitrary multiples of m, 
then the obtained Ai is a wavelet matrix as well. However, if we make these movements 
chaotically, then in general deg(^i) does not coincide with ord(^i) and Ai becomes 
singular. 

As A and UA, where U is nonsingular, have the same rank, order and degree, we 
will additionally assume that A G WM. (m, N, N, F) satisfies (EJ and the subset of 
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such compact wavelet matrices will be denoted by WAio^m, N, N,¥). As has been 
mentioned above, the linear condition does not lose generality in the description 
of WM{m,N, N,¥). 

Matrix polynomials V(z) G WAiolm, 1, 1,F) are called primitive wavelet matrices 
and they always have the form V(z) = Q + Pz, where P and Q are complementary 
(orthogonal) projections on F™' and rank(P) = 1 rank((5) = m — 1 (see ^ Th. 3.1] 
or Lemma 2). 

Every A{z) G WAio{m, N, N,¥) can be uniquely factorized as (see ^ Theorem 
4.4.15] or Theorem 2) 

N N 

(7) A{z) = l[V,{z)=l[{Q,+P,z), V,{z)eWMo{m,l,l,¥), 

i=i j=i 

where no consecutive operators Pj are orthogonal to each other, PjPj^i ^ 0. This 
factorization gives rise to the map (parametrization) 

(8) Fpm-i ^ X ... X FP'"-^ D B < — >WMo{m,N,N,¥), 

N 

where FP™~^ is the projective space F'"/(F\{0}) (the space of one-dimensional sub- 
spaces of F*") which is one-to-one and onto but defined only on the subset B where 
consecutive directions are not orthogonal, i.e. some singular points are excluded from 
the set of parameters. Until now it has been the only known simple way of con- 
structing A G yVAio{m, N, N,¥) for arbitrary A^: choosing nonzero column vectors 
Vj G F™, j = 1,2, . . . , N, such that Vj / "^^+1, taking the corresponding projections 
Pj = Vj{v*Vj)~^v* and primitive wavelet matrices Yj = Im — Pj + PjZ, and constructing 
the product ([7]) that belongs to WMoim, N, N,¥). 

Remark 1. The representation ([^ resembles the factorization for ordinary polyno- 
mials p{z) = Y[k=ii^ ~ ^k) by using their roots. However the following remarkable fact 
should be emphasized: to find the roots of p{z) we have to go, in general, to a larger 
field than the field of coefficients of p{z). For the factorization in IQ, the coefficients 
of each factor Vj belong to the same field ¥. 

Now we turn to our contribution in the study of compact wavelet matrices. This 
approach has been developed during a search for a new matrix spectral factorization 
algorithm [5], and it makes possible to parameterize yVAi{m, N, N, F) directly in terms 
of points in F^"^"^^^. Moreover, compact wavelet matrices can be constructed in a 
more efficient way than by taking the products ([7]). Other advantages and possible 
applications of the proposed method are mentioned in the course of discussions below. 

Recall that we identify A G WM (m, N, N, F) with its polyphase representation 
(j6]). The linear condition we have introduced is ([5]) and such a subclass is denoted 
by yVMlo{m, N, N,¥). We further require that the last row of be not the zero 
vector and denote such a subclass by WA^i(m, N, N, ¥). This is done without loss of 
generality as An ^ and we can interchange the rows of A, if necessary. In particular, 
for any A G WA^(m, A^, A^, F) there exist constant unitary m x m matrices Uq and 
Ui such that UqAUi G >VA^i(m, A^, A^, F) (by right multiplication we assume that 



AUi = {AqUi, AiUi, . . . , AnUi), i.e. the polyphase matrix of AUi is A{z)Ui). Indeed, 
we can interchange the rows of A, if necessary, by multiphcation by Uq and then take 
Ui = [UqA{1)^ . Hence, in what follows, wc parameterize WA^i(m, A*", A*", F). 

Let Pj^^[F] := {J2k=o '^^z^ : Cq, Ci, . . . , c^v G F} be the set of polynomials with coeffi- 
cients from F and PtvI^] '■— { ^k=i ^kZ~^ '■ Ci, C2, . . . , c^v € Fj which can be naturally 
identified with F^+^ and F^, respectively (note that V%^] n T'^fF] = {0} according 
to our notation). Sometimes [F] is omitted because it is clear from the context. 

Let VUi{m, N,¥) be a set of m x m paraunitary matrix functions U{z), 



(9) 

of the special form 



U{z)Uiz) 



(10) lJ{z) = 



^12(2;) 



U2m{z) 
^mm('2') / 



UiAz) e run 



\ Uml{z) U^2{z) 

i ^ i,j ^ TTi, with determinant 1, 

(11) detU{z) = l, 
and with the linear restriction 

(12) U{1) = Im, 

and such that not all polynomials Umj (these are the adjoint polynomials of the entries 
in the last row) are zero at the origin, 



(13) 



E 



mj 



> 0. 



Then for each A{z) e WMi{m, N, N, F), we have U{z) = diag[l, . . . , 1, z-^]A{z) e 
VUi{m, N,¥) (the last row is multiphed by z~-^), and conversely. Thus there is a 
simple one-to-one correspondence 

(14) yVMi(m,N,N,¥) ^VUi(m,N,¥), 

and we parameterize the latter class using the following 

Theorem 1. Let N > 1. For any Laurent matrix polynomial F[z) of the form 



(15) 



F{z)^ 



( 1 









1 









1 











\Cl(^) C,2{Z) C^{Z) 






1 

Cm-i(^) ly 



Q{z) e P^[F], 
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there exists a unique 
(16) 

such that 
(17) 



U{z) e VUi{m,N,¥) 



F{z)U{z) e P+(m X m) 



Conversely, for each U{z) satisfying ff76|) . there exists a unique matrix function F{z) 
of the form ff73]) such that ( [77[ ) holds. 

V^{m X m) stands for the set of m x m matrix polynomials with entries from V~^. 



Remark 2. The relation ( [77[ ) written in the equivalent form U{z) = F~^{z)M{z), 
where M{z) G V^{m x m) and each Q is replaced by —Q in (1731) to get F~^{z), means 
that we have the Wiener-Hopf factorization of the unitary matrix function U{t) = 
U{z)\z=t defined on T, since M{z) and F^^{z) are analytic regular matrix functions 
inside and outside {including infinity) T, respectively. Although we heavily exploited 
this idea in our previous works on the factorization of matrix functions \T\, [S], from 
which this paper has stemmed out, we tried to reduce to a minimum the application 
of Wiener-Hopf factorization theory in the presented proofs and discussions. This 
allowed us to transfer the obtained results to arbitrary subfields F C C which are 
invariant under conjugate operation, and thus extend the area of their applications. 

The proof of Theorem 1 makes it possible to explicitly construct the corresponding 
U{z) for a given F{z), and vice versa. Namely, let the functions in f[T^ be 



N 



likZ 



l,2,...,m- 1, 



k=l 



and let 6, be the (A^ + 1) x (A^ + 1) Hankel-like matrix 



(19) 

and 
(20) 



e,: 



( 


7ii 


1-a ■ ■ 


■ li,N-l 








7i3 ■ ■ 


liN 









7i4 ■ ■ 








\CfiN 





•• 





J 



z = l,2, 



, m - 1, 



m— 1 



i=l 



Assume also that Bi is the first column of 6j, i = l,2,...,m — 1, and = 
(1,0,..., 0)"^. Let Xj = {ajo, aji, . . . , ctj^)'^ G F^+^ be the solution of the follow- 
ing linear system of algebraic equations 

(21) A-X = Bj, j = l,2,...,m 

(A is positive definite and has a displacement structure of rank m, therefore 0(mA^^) 
operations are required for its solution instead of traditional 0(A^^); see P, Appendix]), 
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and let 

N 

(22) Vmj{z) = ^ajkZ-'', j = l,2,...,m, 

k=0 

(23) Vij{z) = [Ci{z)vnij{z)]^ - 6ij, z = 1,2, . . . ,m - 1, 
where [ ■ ]^ stands for the projection operator, 

(24) U{z) = V{z){V{l))-\ 

where V{z) is the m x m Laurent polynomial matrix V{z) = {vij{z))™'j^^, will be the 
desired matrix polynomial f|T6|) as proved in Sect. 4. 

For a given U{z), the corresponding -F(2) can be also explicitly computed as follows 
(see Sect. 4 for the proof): if Mmj(O) ^ 0, then 

i = 1, 2, . . . , m — 1, 

where [ • ]^ stands for the projection operator: [Y^'^=-n '^k^''] = J2k=-N '^kz'^ , and 
under — K-^ its formal series expansion in a neighborhood of is assumed. Note that 
we need to take only the first + 1 coefficients in this expansion in order to compute 

In consequence, as the set of matrix polynomials of the form ( ITSj) can be easily 
identified with F^"^"^-*^, we have the following diagram of one-to-one and onto maps 

(26) WM 1 (m, AT, N, F) < — > VUi (m, N, F) i — > F^ x F^ x . . . x F^^ , 

m—l 

i.e. a complete parametrization of nonsingular compact wavelet matrices, which can 
be effectively reahzed. 

The paper is organized as follows. In the next section we provide all notation and 
definitions used throughout the paper. Most of them have already been introduced 
in the current section but we collect them together for the convenience of reference. 
In Sect. 3 we reprove the factorization (JTj) and related facts as it seems that our 
approach is somewhat simpler than the traditional one. The main result of our paper. 
Theorem 1, is proved in Sect. 4. The essential part of this theorem, the existence of 
the paraunitary matrix function ( !T6|) (without showing the property f|T3|) ). has already 
been published (see [5l Th. 1]) for the case F = C|j It is sufficient to note that the 
proof goes through without any change if we replace C by F. However, for readers' 
convenience, we included the simplified version of this proof in the present paper. (A 
minor inaccuracy of the discussion in [5] between the formulas (50) and (51) has also 
been corrected.) 

In the remaining two sections we consider some possible applications of our method. 
In particular, in Sect 5, we solve the following important problem: given the first row 
of a wavelet matrix A (sometimes called the scaling vector or low-pass filter), how 



E 



-N 



CkZ" 



EN k 
k=0 ^kZ 



Then 



(25) 



Uij{z) 



This constructive proof is a core of a new matrix spectral factorization method published in [5] . 
This method is currently patent pending. 
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to find the remaining rows of A (called tlie wavelet vectors or high-pass filters). Tlie 
solution of this problem is well known and there exists an appropriate algorithm for 
construction of wavelet vectors (see e.g. [3 p. 64]). However this algorithm requires 
too many matrix multiplications which might cause the round-off problems for large 
m and N. We propose another algorithm for solution of this problem which uses our 
approach. A complete comparison of both algorithms based on numerical simulations 
will be provided soon. 

Sect. 6 contains some discussions about other possible applications of the proposed 
parametrization of compact wavelet matrices. 

2. Notation and Preliminaries 

The fields of rational, real, and complex numbers are denoted by Q, M, and C, 
respectively, and F stands for a subfield of C which is invariant under the complex 
conjugation, a G F ^ a G F. 

T := G C : = 1}, T+:= {zeC: \z\ < 1}, and T_ := G C : 1^1 > l}U{cx)}. 

Let ei, 62, . . . , be the standard basis of F™, e.g., ei = (1, 0, . . . , 0). The usual 
scalar product and norm in the space F™ are denoted by ( ■, ■)m and \\-\\m, respectively 
(they are independent of F as it is a subfield of F). 

Let F™"^"^ be the set of m x n matrices with entries from a field F. For M G F™"^", let 
M* = M G F"^*^ be the Hermitian conjugate of M. diag[ai, 02, ... , cim] is an m x m 
diagonal matrix with corresponding entries on the diagonal and = diag[l, !,...,!] 
is the m X m unit matrix. U G F™-^"^ is called unitary, U G Wm[F], if UU* = U*U = Im- 

If M G F*"^"* and a is an entry of M, then cof(a) G F is the cofactor of a and 
Cof(M) G F'"^™ is the cofactor matrix, so that M'^ = ^(Cof(M))^ if M is 
nonsingular. The same notation will be used for matrix functions. 

A matrix P G F"^^"^ (considered as a linear map from F™ to F™) is called the 
(orthogonal) projection if it is self- adjoint, P = P*, and idempotent, P^ = P, however 
we will reduce this definition to a single formula: P = PP* . Note that if P G F*"^"* 
is a projection of rank r < m, then there exists U G C'"^'" with orthonormal columns, 
U*U = Ir, such that 

(27) P = UU* 

(indeed, the factorization fl27p with U G C™^^ of full rank r exists since P is non- 
negative definite, and UU*UU* = UU* U*UU*UU*U = U*UU*U =^ {U*U)-^ x 
{U*U){U*U){U*U)(U*U)-^ = {U*U)-\U*U){U*U){U*U)-^ U*U = Ir). 

If P and Q from F™^™ are projections and P + Q = Im, then they are called com- 
plementary to each other. Obviously, PQ* = Q*P = for complementary projections 
(as P{Im — P)* = P — PP* = 0), and if P is a projection, then — P is also the 
projection complementary to P. 

V\¥] denotes the set of Laurent polynomials with coefficients from a field F, and 
Vn[¥] := {Ef=-7vCfc^'' : Ck e ¥,k = -N,...,N}. If we write just V, the field of 
coefficients will be clear from the context. gV is the set of polynomials (with non- 
negative powers of z, Ckz'' G V^) and C V is the set of Laurent polynomials 
with negative powers of z, X]fc=i '^^-^ ^ ^ • emphasize that, according to our 
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notation, constant functions belong only to so that (1 V = {0}. Let also 

P± = P± n Vn. 

V{m X n) denotes the set of m x n (polynomial) matrices with entries from V, 
and the sets V^{m x n), P^[F](m x n), etc. are defined similarly. The elements of 
these sets, P{z) = (^pij{z)), are called (polynomial) matrix functions. When n = 1 we 
have the vector functions and such classes are denoted by V{m) instead of P(m x 1). 
When we speak about the continuous maps between these sets, we mean that they are 
equipped with a usual topology. 

For p{z) = Y.k=-N'^kz'' e V, let p{z) = Y.k=-N'^kZ~^ and for P{z) = [pij{z)] E 
V{m X n) let P{z) = [pijiz)]^ E V{n x m). Note that P{z) = (P(^))* when z eT. 
Thus usual relations for adjoint matrices like Pi + P2{z) = Pi{z) + P2{z), PiP2{z) = 
P2{z)Pi{z), etc. hold. _ 

We employ also the additional notation of the sets := {p{z) : p{z) E V^[¥]}, 

which is an extension of V^[¥], and V^[¥]{m) := P+[F] x . . . x P+[F] x V^[¥]. 

" V ' 

m 

A polynomial matrix (jH]), where Ak E F"^^"^, k = 1,2, . . . , N, is called paraunitary, 
we write A{z) E VU{m,N,¥), if (i]) holds, which is equivalent to ES ^fc^fc+z = 
I = 0,1,..., N. This means that A = {Aq, A^, As, ■ ■ ■ , A^) E F™>^™(^+i) is 
a wavelet matrix and we do not make distinction between A and its its polyphase 
representation ([6]). If 7^ in the representation ([6]), then we say that the wavelet 
matrix A = A{z) has rank m and order A^, we write A{z) E WAi{m, N,¥). The 
degree of a wavelet matrix A{z) is the order of a monomial det A{z), i.e. det A{z) = 
c-z'^'^^^-^\ The set of wavelet matrices of rank m, order N, and degree d will be denoted 
by yVJ^{m, N,d,¥). In addition, WAio{m, N,d,¥) is the subset of wavelet matrices 
which satisfy (jS]), and A{z) E WMi{m, N, d, F) C yVMo{m, N, d, F) if the last row of 
An differs from the zero row vector. 

Recall also that VUi (m, A^, F) denotes the set of specific paraunitary matrix func- 
tions which satisfy (P])-( TT^ . 

For power series f{z) = Efcl-ooCfe^'' and > 1, let [fiz)]-, [fiz)] + , [/{z)]]^, and 

[/(z)]^, denote respectively Y.kl_^CkZ^, Y.k'=oCkz'', Y.kl-NCkZ^^ and ^^=0^^^:'' and 
the corresponding functions are assumed under these expressions if the convergence 
domains of these power series are clear from the context (the dependence of coefficients 
Cfc on the function / is sometimes expressed by Cfc{/} = c^). Obviously [/]^ = / — [/] + 
under these notation. 

A matrix function S{z) E P[C](m x m) is called positive definite if S{z) is such 
{XS{z)X* > for each 7^ X G C^^™) for almost every z E T. The polynomial 
matrix spectral factorization theorem (see e.g. |3]) asserts that every positive definite 
S{z) E Vn{^ X nT-) can be factorized as 

(28) S{z) = S+{z)S+{z), z E C\{0}, 

where E V^{m x m) and det S~^{z) ^ for each z G T+ (consequently, 5"+ (2;) is 
analytic and invertible in T_), and the representation (125]) is unique in a sense that 



if S{z) = Si{z)Si{z) is another spectral factorization of S{z), then S^i^z) = S^{z)U 
for some unitary U G Um- 

3. Wavelet Matrix Factorization Theorem 

The material of this section is mostly well known, however we provide compact 
proofs of the given statements making emphasis on arbitrariness of a field F. 

Lemma 1. Let 

N 

(29) A{z) = Y,Akz\ AfcGF'"^™, A^^O, 

fe=0 

he a wavelet matrix of rank m, order N and degree d, A{z) G VVA^(m, A^, d, ¥). Then 

a) d> N; 

b) d = N if and only z/rank(ylo) = m — 1. 
Proof. We have 

Y^Alz-' = A{z) = A-\z) = -r-^ACof A{z)f = cz-'{Coi A{z)f , 



k=0 



N ^ 

and since Cof A(2;) G V^{mxm) and^Ar ^ 0, the equation ^ Alz'''= ^^(Cof A(2;)) 

implies a). Using the same reasoning, we conclude that rank(ylo) < m — 1 <^=^ 
Cof(Ao) = ^ d> N (note that if Cof A(2) = Ef=oCfc^^ then Cq = Cof(Ao)). 
Hence b) follows as well (recall that rank(y4o) 7^ m since AqA*^ = 0). □ 

The following lemma will be used in the sequel only for d = 1. 

Lemma 2. (cf. [6, Th. 3.1]). Let V{z) = Q + Pz, P,Q e F™^"^, be a matrix 
polynomial. Then it is a wavelet matrix of rank m, order 1 and degree d satisfying 
V{1) = Im, V(2;) G yVAio{m,l,d,¥), if and only if P and Q are complementary 
projections on F™ and rank(P) = d. 

Proof. The first part of the lemma can be proved by observation that {Q + Pz){Q* + 
P*z) = and V(l) = Im^ PQ* = 0, PP* + QQ* = Im, and P + g = ^ P 
and Q are complementary projections. 

Now the property of P to be a projection is equivalent to the existence oi U E C™^'^, 
rank([/) = rank(P) = d, with orthonormal columns, U*U = Id, such that f l27p holds. 
Let Uq G C™^™ be any unitary matrix which completes the columns of U. Then 

UoV{z)Uo = Uo{I - UU* + UU*z)Uo = diagjz, z, 

d 

whose determinant is z"^. Hence the lemma is proved. □ 

Theorem 2. (see [3 Th. 4.4.15]). Let A{z) be a wavelet matrix of rank m and of 
order and degree N satisfying the linear condition A{z) G yVA4o{m, N, N,¥) . 
Then there exists a unique factorization . 
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Proof. By virtue of Lemma 1, rank(Ao) = m — 1. Let a nonzero column vector v G F™' 
be orthogonal to the columns of Aq, v*Aq = 0, and Vi(2;) = 1^ — P + Pz, where 
P = v{v*v)~^v* is the projection. Then Vi(z) G yVAio{m, 1, 1,F) and yi{z) divides 

A{z) on the left (the quotient B{z) := V^^{z)A{z) = Vi{z)A{z) G V+ (m X m) as 
P*Ao = 0). 

As deg(B) = deg(A)— deg(Vi) = A^— 1 and deg(B) > ord(B) by virtue of Lemma 1, 
we have ord(B) = A^— 1 (obviously, it cannot be less than A^— 1 as Vi(z)B(z) = A{z)). 
Hence B(z) G WJ^o{m, N — 1, A^ — 1,F) and if we continue these divisions, we get 
the factorization ([7j). 

The uniqueness of a divisor Vi(z) (and hence of other factors) follows from the 
fact that if Vi^z) = Q + Pz G WAio{m, 1, 1,F), where P and Q are complementary 
projections and rank(P) = 1 (see Lemma 2), then PAq = as V]"^(z)A(2;) = {Pz~^ + 
Q)A{z) does not have the coefficient at z~^, and since rank(ylo) = m — 1, such a 
projection P is unique. □ 



4. Proof of the Main Result 



We introduce the following system of m conditions which plays a key role in the 
proof of Theorem 1 (where this system comes from is explained in ^ Lemma 5]) . 
Namely, for given Laurent polynomials 



(30) 
let 



(31) 



Qiz)eV^[¥], J = 1,2, 



, m 



(i{z)xm{z) - xi{z) G 

C2iz)Xmiz) - X2{z) G 



< : 



(m^l{z)Xm{z) - X^i{z) G P+, 
^(l{z)Xi{z) + (2{Z)X2{Z) + . . . + (m-l{z)Xm-l{z) + X^{z) G P+. 

We say that a vector function 



1,2, 



(32) n{z) = {ui{z),U2iz),...,Um-i{z),Um{z)) , Ui{z) e V^[¥], i 

(we emphasize that the first m — 1 entries of fl5^ belong to 'P^[F] and the last en- 
try belongs to "P^yf^]) ^ solution of (131]) if and only if all the conditions in (13T]) 
are satisfied whenever Xi{z) = Ui{z), i = l,2,...,m (it is assumed that x^{z) = 
u^i{z)). Observe that the set of solutions of (!3T!) is a linear subspace of V^\¥]{m) := 

P^[F] X ... X Pj^[F] X 'Py^[F]. We will see that actually this subspace is always m 
V ' 

m 

dimensional. 

We make essential use of the following 

Lemma 3. Let i\30h hold, and let 



(33) 



u{z) = {ui{z),U2{z), . . . ,Um-i{z),u^{z))'^ G V^[¥]{m) 
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and 

(34) ^r{z) = {v^{z),V2{z),...,Vm-l{z),^{z)f eV^[¥]{m), 

be two {possibly the same) solutions of the system ([31\\ . Then {u{z),\{z))m is the same 
for every z G C\{0}, i.e. 



m—l 



(35) ^ Ui{z)Vi{z) + Urn{z)Vm{z) = Const . 



i=l 



Proof. Substituting Xi = Vi in the first m — l conditions and Xj = Ui in the last 
condition of fl?T]) . and then multiplying the functions in the first m — l conditions by 
Ui and the function in the last condition by Vm, we get 



C2VmU2 - V2U2 e 



^ClUlVm + C2U2Vm + • • • + Cm^lUm^iVm + U^Vm G 

Subtracting the first m — l functions from the last function in the latter system, we 
get 

m—l 

(36) ^ Ui{z)vi{z) + u;ri{z)v„,{z) e V^. 

i=l 

We can interchange the roles of u and v in the above discussion to get in a similar 
manner that 

m—l 

(37) ^ Vi{z)Ui{z) + V:;n{z)Ujn{z) E . 

1=1 

It follows from the relations fl36|) and fl37|) that the function in fl35|) and its adjoint 
belong to V~^, which implies f l35|) . □ 

Corollary 1. //(E3) ^^'^ (O ^'^^ ^'"'^ solutions of the system (EI]); and 

(38) u(l)=v(l), 

(39) u{z) = v{z) for each z G C\{0}. 

In particular, if u(l) = G F™, then u{z) = /or each z G C\{0} smce the trivial 
vector function v{z) = is always a solution of the system i\31\i . 



Proof. Since u{z) — v(z) is also a solution of fl3T|) . it follows from the lemma that 
||u(z) - v{z)\\m = ||u(l) - v(l)||^ = for each z G C\{0}. Hence ([39]) holds. □ 

Remark 3. One can see that Corollary 1 remains valid if we take any fixed point 
zq ^ in the role of 1 in the equation i\38^ . 
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Lemma 4. Let ( Tf^ hold, and let vector functions 

(40) ^rj{z) = {vi,{z),V2,{z),...,^j;Zj{z)f eV^[¥]{m), j = l,2,...,m, 

6e an^/ m solutions of the system {\31\) . Then the determinant of the m x m Laurent 
polynomial matrix 

(41) V(z) = (vi(z),V2(z),...v^(z)) 
is constant 

(42) det V{z) = Const ., ze C\{0}. 

Proof. Since the vector polynomials (l40l) satisfy the last condition of the system (l3T!l . 
we have 

(43) F{z)V{z) e P+(m x m), 
where F{z) is defined by f|T5|) . Consequently, as detF(z) = 1, 

(44) detl^(t)eP+. 

Since the vector polynomials (I40p satisfy the first m — 1 conditions of the system (1311) . 
we have 

(f)ij{z) := Ci{z)vmj{z) - vtj{z) e V^, 
l<j<m, l<i<m. Direct computations show that 

V{z)F-\z) = $(2) G P+(m X m), 

where F~^(2;) is obtained from F{z) by replacing each Q by — and 



( -(plliz) -021 (2:) ■■■ -(t)m-l,l{z) Vynl{z)\ 
-012 (z) -022(2;) ■■■ -0^-1,2(2;) Vyn2{z) 



\-0ml(2) -0ml (2;) ■■■ -0m-l,m(2;) Vm,m{z) / 

Consequently det i/(z) = det ^{z) E and together with (jH]) this gives (j42]). □ 

Remark 4. Observe that the relation ( |^5D /lo/ds under the hypothesis of the lemma. 

Now we construct explicitly m independent solutions (140|) of the system (!3T|) . We 
seek for a nontrivial solution 

(45) x(z) = (xi(2),X2(2),...,5;;(z))'^GP®[F](m) 
of (pTI) . where 

N 

(46) Xi(z) = ^aifc^^, i = l,2, ...,m, 

fc=0 

and the coefficients ajfc are to be determined. 

Equating the coefficients of z~^, k = 0,1,2, N of the Laurent polynomials in 
the system fl3ip to 0, except the 0th coefficient of the jth function which we equate to 
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1, we get the following system of algebraic equations in the block matrix form which 
we denote by Sj: 

©l-^m —-^1 = 0, 
— X2 = 0, 

(47) Sj:= {QjX^-X]=l, 



GiXi + Q2X2 + . . . + &m-lXm-l + Xm — , 

where 6^ is defined from the equations ^ and = (0,0,..., 0)^ G F^+\ 

1 = (1, 0, 0, ... , 0)"^ G F^"*"^, and the column vectors 

(48) Xi = {aio,aii,...,aiN)^, z = l,2,...,m, 

(see (H6|) ) are unknown. 

Remark 5. We emphasize that if (Xi, X2, . . . , Xm) defined by i \48^ is the solution of 
the system (^7), then the vector function ( |^5D defined by ( [7^ will be the solution of 



the system i[3lh . 

It is easy to show that the system fll7|) Sj has a solution for each j = 1,2, ... ,m. 
Indeed, determining Xi,i = l,2,...,m — 1, from the first m — 1 equations in ( H71) . 



(49) Xi = ei-X^-6ijl, 2 = l,2,...,m-l, 

and then substituting into the last equation of (1471) . we get 



Ol • Ol ■ Xm + O2 ■ ©2 ■ Xfn + . . . + Om-l ' ©m-l ' X^ + -^m — ' 1 

(we assume that 9^ = ^at+i, i-e. the right-hand side is equal to 1 when j = m) or, 
equivalently. 



(50) (Gi ■ Gi + 62 ■ 62 + . . . + Qm-i ■ <S>m-i + In+i) ■ Xm = ■ 1. 

This is the same system as ( |2T!) (see [20]). Since each 9j is symmetric, 9i = Of, 
we have that Oj0j = OjO*, i = 1,2, ...,m — 1, are non-negative definite and the 
coefficient matrix (1201) of the system (l50l) (which is the same for each j = 1,2, ... ,m) 
is positive definite (with all eigenvalues larger than or equal to 1). Consequently, A is 
nonsingular, det A > 1, and the system fl50l) has a unique solution for each j. 

Finding the vector Xm G F^"*"^ from ( l50l) and then determining Xi,X2, ■ . . ,Xm-i 
from ( l49l) . we get the unique solution of Sj. To indicate its dependence on j, we denote 
the solution of by (X^, X^, X^.^, X^), 

(51) XI := {al^,al^,...,alj^f, i = 1,2, . . . ,m, 
so that the vector functions ( l40l) . where 

TV 



k=0 
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are m solutions of the system ( 1311) (see Remark 5). These vector functions vi{z), V2(-2), 
. . ., Vm{z) are linearly independent since the linear transformation L : 'P^[F](m) — )■ F™ 
which maps (xi(z), X2{z) , . . . , x^{z)) into the 0th coefficients of the functions (or their 
adjoints) standing on the left-hand side of the system ( 13T|) . viz. into {co{Ciiz)x^(z) — 

xi{z)}, Co{C2{z)x^{z) -X2{z)}, . . ., Co{Ci{z)xi{z) + (2(2) X 2(2) + ■ ■ . + x^{z)}) , trans- 
forms m vector functions vi(z), V2(z), . . ., "Vm{z) into linearly independent standard 
bases of F'". Namely, L(vj(z)) = {6ii, 6i2, . . . , Sim), i = 1,2, ... ,m, because of (H71) . 
Consequently V{1) is nonsingular (see fHTj) ) since if w ■ V{1) = G F"* for some 
7^ w = {wi,W2, . . . ,Wm) £ IF™, then w ■ V{z) = J2T=i '^k^ki.z) = by virtue of 
Corollary 1 of Lemma 3, which contradicts the independence of ^i{z), \2{z),. ■ .,Vm{z). 
Let U{z) be defined by fl24|) . Then it has the form (|T0|) and its column vectors ui{z), 

U2{z),. . .,nm{z), 

(52) U{z) = {m{z),U2{z),...,Um{z)), 
are m solutions of (13T!) satisfying 

(53) Ui(l) = ei, i = 1,2, . . . ,m, 

(since f|T2l) holds because of fl24l) ). Therefore, the matrix function U{z) is paraunitary 
since U{z)U{z) is a constant matrix, by virtue of Lemma 3, and this constant is /„, 
since (fT2|) holds. By virtue of Lemma 4, the determinant of U{z) is also a constant 
which is equal to 1 since detM(l) = 1. Hence the relations ©-([I2D are valid for 
(15^ . Observe also that (ITTIl holds because of (see Remark 4) and (^^. Hence 
it remains to show only the relation ( fT3l) in order to complete the proof of the first 
part of Theorem 1 . Meanwhile we are ready to formulate and prove a corollary of the 
above discussion which will be used in the next section. 

Corollary 2. Let {\3(^ hold, let F{z) be the matrix function defined by {\15\] . and let 
ff?^ be the corresponding {according to Theorem 1) U{z) e VUi{m, N,¥) . 

J/b = (&i, b2, . . . , fem) ^ then u^y^z) is a solution of the system (f371) satisfying 

Ub(l) = b, 

if and only if 

m 

Ui,{z) =J2biMz), zeC\{0}. 

i=l 

Proof. Since ui{z), U2{z), . . . , Um{z) (see ([S2D) are solutions of the system (IHT]) (recall 
that the set of its solutions is a linear subspace) and fl53p holds, the first part of the 
corollary is clear. The second part follows from the first part and from Corollary 1. □ 

In order to prove f|T3|) . let 

(54) ^{z) = {i,,,{z))l^^:=Fiz)Uiz), 

which, as has already been mentioned, belongs to V^{m x m) (see (fT71) ). Since the 
determinants of F{z) and U{z) are 1 for each z G C\{0}, we can conclude that 

(55) det "^{z) = 1 for each z e C. 
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Because of the structure of the matrix ffTSl). we have 



Uij[Z) 



1 < i < m, 1 < j < m. 



Since U{z) = U ^{z) and ( ITTi) holds, we have 

Urnj{z) = Cof Umj{z) = Cof 1pmj{z), I < ] < TTl. 

Thus, by virtue of fl55|) . 

m m 
1 = det*(z) = ^^/'mil^) cof = y^^iJmj{z)Umj{z) 

for each 2; G C including 0. Consequently ffT^ holds. 

Thus, the constructed matrix function U [z) has all the desired properties, i.e. f ll6p 
and f[T7|) hold. Note that a brief way of construction of the entries of U{z) is described 
by the formulas ([T5])-(|2I]). 

The uniqueness of (fT6l) follows from the uniqueness of spectral factorization (see 
Section 2) since F{z)U{z) is the spectral factor of F{z)F{z). 

Let us now show the converse part of Theorem 1. If we have a matrix function 
f|T6|) . then Umj(O) 7^ for some j since f|T3l) holds. Let us determine functions Qi by 
the formula fl2^ . construct the matrix function f[T^ . and show that f[T7|) is valid. For 
this we need only to check that the entries of the last row of '^{z) (see flM|) ) belong to 
[F] . For 1 < n < m, we have 



m— 1 

P 9 ^0(^)Mm(2) +Ur, 



i=l 



m— 1 



Uij{z)_ 



' m— 1 



m— 1 r 



(56) 



Umji^z) 



Uij{z)Uin{z) + Umj{z)Umn{z) ) - ^ 



Umji^z) 



Uin{z) = 



m— 1 

E 

i=l 



Mi 



which is analytic in a neighborhood of and yields that this function belongs to V^. 
Thus (dZD holds. 

Let us now show the uniqueness of the desired F{z). 

Recall that F~^{z) can be obtained by replacing each Q by —Q in f[T^ . Hence 

F-'^{z) e P+[F](m X m). Since ^(^) G P+(m x m) (see (EID and ([IT])) and holds, 
we have \l/^^(z) G P^(m x m). 

If Fi(^) has the same form f|T5l) with the last row [Ci(2;), C2(^)5 • • • 5 C'm-ii^)^ 1] 
^^(^) := Fi{z)U{z) G P+(m x m), then U{z) = F-^{z)-^{z) = F^^{z)^i{z), which 
yields 

V+{m xm) 3 Fi{z)F-^{z) = ^i(z)$"^(2) G P+(m x m). 

Thus Fi(z)F~^(z) is a constant matrix while, on the other hand, this product has the 
form ( [T5|) with the last row [([ — (i, (2 — (2, ■ ■ ■ , C'm-i ~ Cm-i, 1], which implies that 
C = Ci, ^ = 1,2, . . .,m - L 
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5. Wavelet Matrices with Given First Rows 

If we know the first row of the wavelet matrix ([71), then one can uniquely determine 
each Pi and Qi in (JTj) and then recover A.{z) by the formula ([7]), i.e. find the remaining 
m — 1 rows of A{z). The algorithm for this procedure is well known (see e.g. [71 
Th. 4.4.17]). We describe a new method of reconstruction of A(z) based on the 
proposed parametrization of compact wavelet matrices. This description is included 
in a constructive proof of the following 



Theorem 3. Let v = {vi,V2, ■ ■ ■ , Vm) G F™ be a vector of the unit norm and V G F™^™ 
be a unitary matrix with the first row v. If Ai{z) = (^au{z) , ai2{z) , . . . , aim{z)) E 
P^(l X m) is a row function of order N (^^i \ciNj\ > 0, where aij{z) = ^ 
j = 1, 2, . . . , m) such that Ai{z)Ai{z) = 1^, z G C\{0} and 



N k 
k=0 ^kjZ , 



(57) 



Ai(i; 



then there exists a unique wavelet matrix of rank m and of order and degree N, A{z) G 
yVM{m, N, N,¥), with A(l) = V which has the first row Ai{z). 



Remark 6. For a given v G F™ with \\v\\m = I , it is in general hard to decide whether 
there exists a unitary matrix V G F"*^™ with the first row v {because of the arbitrariness 
of the field F which not always contains ^/a for a positive a G F) . However the theorem 
can be always used for v = ei and V = Im- 

Proof. For notational convenience, we will find A'^{z) instead of A{z). Thus we assume 
that the v = (t>i, V2, ■ ■ ■ , is a column vector and Ai{z) = (011(2), . . . , aim{z))'^ G 

V^ijn X 1), which satisfies Ai{z)Ai{z) = 1^ and (1571) . Furthermore, without loss of 
generality (we interchange the rows if necessary), we assume that a^rn 7^ where we 
recall aim{z) = Y!k=o(^kmZ^. 

Consider now Ui(z) = {uii{z),U2i{z), . . . ,u^i{z)) G V^{m), where Uji{z) = 



can define Q, 



5 

i = 



m 



z ^a,im[z) 



1,2, 



1, and Umi{z) 
, m — 1 , by a formula like (125 



Then Umi(O) 7^ 0, and we 



O(^) 



" Uii{z) ' 




' Uii{z) ' 


Uml{.z)_ 




Uml{z)_ 



l,2,...,m- 1, 



N 



and construct the matrix function F{z) defined by ( lT5l) . Consider now the corre- 
sponding (according to Theorem 1) U{z) G VUi{m, N, F) whose columns ui{z), U2{z), 
. . .,Um{z), as we remember, are solutions of the system ( 13T1) and also satisfy (1531) . We 
can make sure by direct computations like dSS]) that Ui{z) satisfies the last condition 
of (1311) . and taking into account the relations 



Ciiz)Umliz)-Uii{z) 



Uil{z) 
Uml{z) 



Uiljz) 
.Uml{z) 



Umliz)-Uii{z) 



Uil 
Uml 



Uml G 
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z = l,2,...,m — 1, we see that Ui{z) is a solution of the system ( 1311) . Consequently, 
according to Corollary 2, 

m 

Ui{z) =U{z) ■v = ^^ViUiiz), 

i=l 

(as this equation holds for z = 1) and it is the first column of U (z) ■ V. This matrix 
function is paraunitary as well (since V is unitary) and satisfies the condition U{1)V = 
V. Hence 

A(z) = diag[l,l,...,l,;2^]f/(z)l^ 

will be the desired matrix. 

The uniqueness of A{z) is also valid since F{z) described in the above proof is the 
same as the matrix function corresponding to A{z) ■ V^^ G WA1i(m, N, N, F) by the 
to one-to-one map (HM and Theorem 1, and such F{z) is unique. □ 

6. Other Possible Applications 

In the end, we would like to discuss briefly some possible applications of the proposed 
method of compact wavelet matrices construction. It should be emphasized that the 
intent of this section is largely motivational and we consider three separate topics. 

A. Rational approximations to compact wavelet matrices. Most of compact 
wavelet matrices used in practice, for example Daubechies wavelet matrices, obey 
certain additional restrictions. The proofs of the existence of such matrices and the 
ways of their construction are highly non-linear and thus the obtained coefficients are 
in general irrational. In actual calculations on a digital computer, these coefficients 
should be quantized and hence approximated by rational numbers. It may happen 
during this approximation that the basic property of the wavelet matrices ([2]) will not 
be preserved exactly. Using the proposed parametrization and taking Q in the role 
of F, we can provide an approximation of any compact wavelet matrix A by A' with 
rational coefficients preserving exactly the shifted orthogonality condition ([2]). Indeed, 
it is just sufficient to note in the proof of Theorem 1 that the maps described in fl26|) are 
continuous. Hence, if we find the point in M('"~i)-'^ corresponding to A, approximate 
it by rational coordinates and go back into the space of wavelet matrices, then we will 
get the desired A' (in a similar manner, Vaidyanathan [8J used the parametrization ([8]) 
for the same purposes). Such various approximations for Daubechies wavelet matrices 
of different genus are explicitly constructed in [1]. 

B. Cryptography. As said in Remark 1, the factorization ([7]) much resembles the 
factorization for ordinary polynomials into linear terms. On the other hand, an efficient 
way of construction of paraunitary matrix polynomials associated to the given points in 
]p(m-i)7v expressed by the diagram fl26l) . which helps to handle such matrix polynomials 
easily, can be compared to natural parametrization of ordinary polynomials by to their 
coefficients. An application of polynomial factorization theory in the cryptography 
is widely known. Having the above similarities (and the advantages mentioned in 
Remark 1), one can also expect certain applications of the developed theory in the 
cryptography. 
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C. Selection of best wavelet matrices. Which wavelet matrix is most suitable 
to apply in a given practical situation represents frequently a certain optimization 
problem (which is sometimes very hard to solve) or should be obtained empirically 
by computer simulations. Having a quick access to the complete bank of compact 
wavelet matrices due to the parametrization ( 126|) gives an opportunity to choose the 
best possible wavelet matrix by nearly complete screening. 

References 

[1] L. Ephremidze, A. Gamkrelidze, and E. Lagvilava, "An approximation of Daubechies wavelet 
matrices by perfect reconstruction filter banks with rational coefficients" , (to appear in Adv. 
Comput. Math.). 

[2] L. Ephremidze, G. Janashia and E. Lagvilava, "On the factorization of unitary matrix-functions" , 
Proc. A. Razmadze Math. Inst, vol. 116, pp. 101-106, 1998. 

[3] , "A simple proof of matrix-valued Fejer-Riesz theorem", J. Fourier Anal. AppL, vol. 15, 

no. 1, pp. 124-127, 2009, (DOI: 10.1007/s00041-008-9051-z). 

[4] , "On approximate spectral factorization of matrix functions" , to be published in J. Fourier 

Anal. AppL, (DOI: 10.1007/s00041-010-9167-9). 

[5] G. Janashia, E. Lagvilava, and L. Ephremidze "A new method of matrix spectral fac- 
torization", IEEE Trans. Inform. Theory, vol. 57, no. 4, pp. 2318-2326, 2011, (DOI: 
10.1109/TIT.2011.2112233). 

[6] J. Kautsky and R. Turcajova, "Pollen product factorization and construction of higher multiplicity 
wavelets", Lin. Algebra AppL, vol. 222, pp. 241-260, 1995. 

[7] H. L. Resnikoff and R. O. WeUs, Wavelet Analysis, Springer- Verlag, 1998. 

[8] P. P. Vaidyanathan, Multirate Systems and Filter Banks, Prentice Hall, New Jersey, 1993. 

Authors' Addresses: 

L. Ephremidze, E. Lagvilava 
A. Razmadze mathematical Institute 

1. Javakhishvili State University 

2, University Street, Tbilisi 0143, Georgia 
E-mails: lephremi@umd.edu; edem@rmi.ge 



